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ABSTRACT 

This paper presents two ease studies where a maeroseopie model-based approaeh for traffie state 
estimation, whieh we have reeently developed, is employed and tested. The estimation method¬ 
ology is developed for a “mixed” traffie seenario, where traffie is eomposed of both ordinary and 
eonneeted vehieles. Only average speed measurements, whieh may be obtained from eonneeted 
vehieles reports, and a minimum number (suffieient to guarantee observability) of spot sensor- 
based total flow measurements are utilised. In the first ease study, we use NGSIM mieroseopie 
data in order to test the eapability of estimating the traffie state on the basis of aggregated in¬ 
formation retrieved from moving vehieles and eonsidering various penetration rates of eonneeted 
vehieles. In the seeond ease study, a longer highway stretch with internal congestion is utilised, 
in order to test the capability of the proposed estimation scheme to produce appropriate estimates 
for varying traffic conditions on long stretches. In both cases the performances are satisfactory, 
and the obtained results demonstrate the effectiveness of the methodology, both in qualitative and 
quantitative terms. 


Keywords: traffic state estimation, connected vehicles, Kalman filter, real-data testing 



Roncoli, Bekiaris-Liberis, Papageorgiou 


3 


INTRODUCTION 

Traffic state estimation is a task of crucial importance for the development and application of traf¬ 
fic management and control strategies. Essentially, it consists in estimating the traffic variables 
of a highway network at sufficient spatial resolution in real-time based on limited amount of in¬ 
formation available from real-time traffic measurements. In conventional traffic, the necessary 
measurements are provided by spot sensors (based on a variety of possible technologies), which 
are placed at specific highway locations. In case the sensor density is sufficiently high (e.g., ev¬ 
ery 500 m), the collected measurements are usually sufficient for traffic surveillance and control; 
otherwise, appropriate estimation schemes need to be employed in order to produce traffic state 
estimates at the required space resolution; see, for instance, (i, 2, 3, 4) among many other works 
addressing the problem of highway traffic estimation employing conventional detector data. 

In the last two decades there has been a significant and increasing interdisciplinary effort by the 
automotive industry as well as by numerous research institutions around the world to plan, develop, 
test and start deploying a variety of Vehicle Automation and Communication Systems (VACS) 
that are expected to revolutionise the features and capabilities of individual vehicles within the 
next decades. A wide description of VACS, focusing particularly on their impact on highway 
traffic efficiency, may be found in (5). To achieve related traffic flow efficiency improvements on 
highways, it is of paramount importance to develop novel methodologies for modelling, estimation 
and control of traffic in presence of VACS. Several papers are providing interesting results related 
to traffic control in presence of VACS, employing either microscopic or macroscopic approaches; 
see, for example {6, 7, 8, 9, 10,11,12,13,14,15,16,17). 

Some VACS are characterised by different types of communication devices, providing the possi¬ 
bility of V2V (vehicle-to vehicle) or V2I (vehicle-to-infrastructure) communication; specifically, 
“connected” vehicles can report their position, speed, and other relevant information, i.e., they 
can act as mobile sensors; existing examples include mobile phones and network-connected GPS 
(Global Positioning System) devices. From a traffic estimation viewpoint, this may allow for a 
sensible improvement of the achievable estimation accuracy, as well as significant reduction of the 
necessary number of spot sensors; the latter would lead to a sensible reduction of the purchase, 
installation, and maintenance cost of traffic monitoring. 

The exploitation of information retrieved from connected vehicles for traffic estimation purposes 
is considered in numerous research works. Work et al. (18) designed an ensemble Kalman filter for 
a Godunov-discretised Lighthill-Whitham-Richards (LWR) model, estimating the velocity field on 
highways using data obtained from GPS devices; a real-case test is defined based on data retrieved 
from the Mobile Century experiment (79). The work by Treiber et al. (20) is based on data fusion 
techniques, utilising information retrieved both from fixed and mobile sensors; two case studies 
are reported, based on data from the highway A9, near Munich, Germany, and from highway M42, 
near Birmingham, UK. Herrera and Bayen (27) designed a methodology to incorporate mobile 
probe measurements into a modified LWR model, exploiting the so-called Newtonian relaxation; 
the efficiency and performance of the estimation algorithm were tested on highway T880 in Cali¬ 
fornia within the frame of the Mobile Century experiment (79). Yuan et al. (22) formulated a traf¬ 
fic model in a Lagrangian coordinate system (where state variables move with the traffic stream), 
including both Eulerian and Lagrangian measurements; an extended Kalman filter is then used 
to obtain appropriate traffic estimates; and the proposed method is tested on data retrieved from 
highway M42, near Birmingham, UK. Piccoli et al. (23) addressed the impact of both sampling 
frequency and penetration rate on mobile sensing of highway traffic, employing a second-order 
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phase transition model (PTM) and proposing several data fusion schemes to incorporate vehicle 
trajectory data into the PTM; extensive tests using NGSIM data (24) are performed. Seo et al. (25) 
proposed a novel probe vehicle-based estimation method for obtaining volume-related variables 
by assuming that a probe vehicle can measure the spacing to its leading one; remarkably, they do 
not assume the knowledge of any fundamental diagram and perform a real experiment employing 
prototypes of probe vehicles. This method is extended and improved by Seo and Kusakabe (26) 
taking into consideration the conservation law. Bekiaris-Liberis et al. (27) proposed a macroscopic 
model-based approach for the estimation of the total density and total flow; i.e., the density/flow 
of all vehicles (connected and ordinary), for the case of “mixed” traffic utilising speed measure¬ 
ments that may be reported by connected vehicles, and a small number of spot flow measurements; 
this approach removes the requirement of (empirical, hence uncertain) traffic speed modelling, 
such as the fundamental diagram. The performance of the developed estimation schemes was val¬ 
idated through simulation using the well-known second-order traffic flow model METANET (28) 
as ground truth. 

In this paper, we present two case studies using real data retrieved from highways in order to test 
and evaluate the behaviour of the recently developed estimation scheme (27). Specifically, we 
address the problem of estimating the (total) density and flow of vehicles in highway segments 
of arbitrary length in presence of connected vehicles. The two case-studies are built using data 
stemming from different real experiments. The first case, based on the NGSIM microscopic data, 
aims at verifying the capability to estimate the traffic state at high space-resolution (i.e., consider¬ 
ing segments in the order of 50 m in length), using aggregate information retrieved from moving 
vehicles, that is, the speeds of connected vehicles and the total flow measured at the entrance and 
exit of the highway stretch; the proposed estimation scheme is also tested for different penetration 
rates of connected vehicles. In the second case, a longer highway stretch and a longer time horizon 
are considered, using information coming from fixed detectors, namely speed measurements at 
every segment (emulating connected vehicles reports) and a small amount of flow measurements. 
This permits to verify the capability of detecting the onset of a congestion created within the net¬ 
work, which is one of the most crucial and challenging aspects in case this approach is exploited in 
conjunction with a traffic control strategy. In both cases, the estimation results are analysed quali¬ 
tatively as well as quantitatively, employing an appropriate metric to compare the performance of 
the proposed methodology under different scenarios. 

In the following section, we present a linear time-varying model of the density dynamics and a 
Kalman filter for density estimation; then two case studies to evaluate the performance of the 
estimation scheme are described; finally, the paper is concluded highlighting the main findings and 
introducing further challenges for future research. 

TRAFFIC ESTIMATION USING AVERAGE SPEED MEASUREMENTS 
Traffic Density Dynamics 

We consider the following discrete-time equations that describe the dynamics of the total densities 
pi, namely the number of vehicles within segment i divided by the the length of the segment A/, on 
subsequent highway segments (see, e.g., (29)) 

pi(k+l)= pi(k) + ^ {qi-i(k) - qi(k) -f r,(k) - Si(k )), 


( 1 ) 
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where i = 1,... ,A^ is the index of the specific highway segment, N being the number of segments 
on the highway. For all traffic variables, we denote by index sub-z its value at the segment i of 
the highway; qt is the total flow at the exit of segment z; T is the time-discretisation step; and 
A: = 0,1,... is the discrete time index. The variables r/ and Sj denote the inflow and outflow of 
vehicles at on-ramps and off-ramps, respectively, at segment z. Using the known relation 


qi{k) =pi{k)vi{k), 


( 2 ) 


where v/ is the average speed in segment z, we write Equation 1 as 

Pi{k+^) = ^Vi-i(k)p,_i(k) + (^1 - Pi(k) + ^ inik) - Si{k)). (3) 

Note that, in order to guarantee that the discrete-time relation described by Equation 2 is suffi¬ 
ciently accurate, the following inequality must be satisfied 

max^vz(k)<l. (4) 

i,k 1 


Assuming that the unmeasured on-ramp and off-ramp flows are constant (or, effectively, slowly 
varying), the unmeasured ramp flow dynamics may be reflected by a random walk, i.e., 

di{k+\) = di{k)^^f{k), (5) 

where is zero-mean white Gaussian noise and 


0 ,= 


I 




if Hi G Lr 
if rii e Ls 


( 6 ) 


for all z = + Is, with L;. = (zzi,..., zz/J and Ls = {m^+i , • • •, nq+i,} being the sets of seg¬ 

ments, denoted by zz,, which have an on-ramp or an off-ramp, respectively, whose flows are not 
directly measured; and ly and Is are the number of unmeasured flows at on-ramps and off-ramps, 
respectively. We assume that at a segment z there can be either only one on-ramp or only one 
off-ramp, which is typically the case on a highway, and hence, fl = 0. 

Assuming that the average speed of vehicles in a segment, namely v/, is measured (e.g. from 
connected vehicle reports) and defining the state 


-^= {pi,---,pN,0i,...,0q+iy, 


(7) 


the deterministic part of the dynamics of the total density given in Equation 3 and of 0; given in 
Equation 5 can be written in the form of a linear time-varying system as 


x{k+ 1) = A{k)x{k) +Bu{k), 


( 8 ) 
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where 


A(k) = 


B = 

u{k) = 


^ij = 1 

^tiii ~ 

«M,7 = — 1, 

(^ij — Ij 
^ (^ij — 0 , 


if / — 7 = 1 and i>2 
if i = j 

if rii G Lr and j =N + i 
if rii G Ls and j = N + i 
if N < i < Ni and j = i 
otherwise 


bij = ^., if/=landj = l 

bmu = if nii^L, I <mi<N,l <i< N 2 , and j 

bij = 0, otherwise 


Ui = qoik), ifi=l^ 

Ui^i = fffi^ S/fip if ni; ^ L 


/ + 1 


(9) 


( 10 ) 

( 11 ) 


with L = LrUL„Ni=N + lr + ls,N 2 =N-lr- h, A G Be (^2+0. Note that go, 

which is assumed to be measured (for example, via fixed flow detector), denotes the total flow of 
vehicles at the entry of the considered highway stretch and acts as an input to Equation 8; along 
with any directly measured on-ramp and off-ramp flows, r,- and 5,-, i ^ L^, i ^ respectively; while 
Vi, i= are viewed as time-varying parameters of Equation 8. 

We now turn our attention to the measured outputs. We assume availability of the density (or, 
equivalently, of the flow) at the mainstream exit of the highway. If there is exactly one unmeasured 
ramp within the considered highway stretch, then no additional measurements are necessary for 
flow observability. On the other hand, if there are more than one unmeasured ramps within the 
stretch, we need one mainstream flow measurement at any highway segment between every two 
consecutive unmeasured ramps. 

In summary, the measured outputs associated with Equations 8-11 are the density (or, equivalently, 
the flow) at the exit of the considered highway stretch and at a highway segment between every 
two consecutive ramps whose flows are not measured. Therefore, 


y{k)=Cx{k), 


where C G M.^b+h)>k{N+ir+is) defined as 


C 


Cij = 1, for all z = 1,..., 4 -b 4 — 1 and some n* < j < — 1 

Cij = 1, if i = Ir + h and j = N 
Cij = 0, otherwise 


( 12 ) 


(13) 


where L* = jzzj, ..., | is the set L ordered by <. 

Although it is physically intuitive that the system described in Equations 8-13 is observable, when 
the additional mainstream fixed flow sensors are placed at some segment between every two con¬ 
secutive unmeasured ramps, we rigorously proved in (27) that the system is indeed observable 
for certain cases, such as, for example, when a fixed sensor is placed on the mainstream at every 
segment immediately before an unmeasured ramp. 

We summarise below the measurement requirements for the proposed estimation algorithm. 
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Traffic System 



FIGURE 1 The traffic system under consideration and the proposed Kalman filter imple¬ 
mentation. The data used to operate the Kalman filter are speed measurements coming from 
connected vehicles (solid lines) and flow measurements coming from fixed sensors (dashed 
lines). The variable denotes the measurement of quantity w at segment i, which might he 
different than the actual quantity w, due to, for example, the presence of measurement noise. 
A variable represents the value of quantity w of connected vehicles, which may be present, 
at segment i. 


• The average speed of eonneeted vehieles at a segment of the highway, v,-, i = 1,..., A, 
are measured; the assumption being that there is no systematie differenee between the 
average speed of eonneeted vehieles and the average speed of all vehieles in a segment. 

• The total flow of vehieles at the entry and exit of the eonsidered highway streteh, qo and 
qM, respeetively. 

• Either the total flow of vehieles at a ramp, say, r/ or Si, is measured or an additional 
mainstream flow measurement, say qj, at any highway segment between two eonseeutive 
unmeasured ramps, is available. 


Kalman Filter 

We employ a Kalman filter for the estimation of the total density of vehieles in eaeh segment of a 
highway (Figure 1). Defining the estimated state veetor 

x= (pi,...,pyv,0i,...,0z,+4)^, (14) 

the Kalman filter equations are (see e.g., (30)) 
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x{k+ 1) 

= A{k)x{k) +Bu{k) + A{k)K{k) {z{k) —Cx{k)) 

(15) 

m 

= P{k)C^ {CP{k)C^ ARY^ 

(16) 

P{k+l) 

= A{k){l-K{k)C)P{k)A{k)'^ AQ, 

(17) 


where z is a noisy version of the measurement y, R = > 0 and Q = > 0 are tuning param¬ 

eters. Note that, in the ideal case in which there is additive, zero-mean Gaussian white noise in 
Equations 8 and 12, respectively, R and Q represent the (ideally known) covariance matrices of the 
measurement and process noises, respectively (in particular, the qN+m+i elements of Q represent 
the filter’s anticipation for the covariance of (^®). Since the system equations here are relatively 
complex, some tuning of R, Q may be needed for best estimation results. The initial conditions of 
the filter described by Equations 15-17 are chosen as 

xiko) = II (18) 

Piko) = H, (19) 

where /i and H = > 0, which, in the ideal case in which x{kQ) is a Gaussian random variable, 

represent the mean and auto covariance matrix of x{kQ), respectively. The Kalman filter delivers 
estimates of the segment total densities p,-, as indicated at the output of the Kalman filter in Eigure 1. 

EXPERIMENTAL TESTING OF THE PROPOSED ESTIMATION METHODOLOGY 
Case Study 1: NGSIM1-80 Data 

Network and Data Description 

In order to evaluate and illustrate the effectiveness of the proposed estimation scheme, microscopic 
traffic data collected within the Next Generation SIMulation program (24) are utilised. However, 
since these data incorporate non-negligible errors in the position of individual vehicles {31), some 
correction methodologies were proposed to improve their reliability {32, 23). This work utilises 
the data processed by Montanino and Punzo {32, 33), which include the trajectories of all vehicles 
travelling along a stretch in the northbound direction of T80 in Emeryville, California, recorded 
from 4:00 PM. to 4:15 PM on April 13, 2005. The highway is composed of 6 lanes; however, lane 
1 is a so-called HOV (high-occupancy vehicle) lane, characterised by access restricted to a limited 
set of vehicles (only vehicles with a driver and one or more passengers are allowed), causing its 
traffic characteristics to be structurally different compared with the other lanes; therefore this lane 
is excluded from our estimation. The considered stretch (sketched in Eigure 2) is 400 m long and 
an on-ramp is entering the mainstream, where the merge nose is located 175 m after the network 
origin. 

As it can be seen from Eigure 3, a massive congestion is present within the stretch, where con¬ 
gestion waves are coming from downstream and crossing the entire stretch. In order to perform 
macroscopic evaluations, the stretch is divided into N = 8 homogeneous segments of 50 m in length 
(see Eigure 2); the on-ramp is placed within segment 4; while a discrete time step T = 5 s is used. 
Even though this configuration slightly violates the condition of Equation 4 for some time inter¬ 
vals and some segments, the obtained results were found to be very similar to those obtained when 
using lower values for the time step T, which satisfy fully the condition of Equation 4. 
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1 2 3 4 5 6 7 8 


FIGURE 2 A graphical representation of the stretch of the highway 1-80 in Emeryville, 
California, used in Case Study 1. 
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FIGURE 3 The aggregated vehicle speeds extracted from the NGSIM data. Congestion 
shock-waves are spilling hack from downstream and cross the entire highway stretch. 
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All necessary information to perform estimation is extracted from the available trajectory data. In 
particular, the segment speed Vi{k) is computed as the arithmetic average of the speeds of con¬ 
nected vehicles present within the segment at the measurement time instant kT ; whereas total flow 
measurements are computed via two virtual spot detectors; one placed at the network entrance 
(beginning of segment i= 1), providing qo, and another placed at the network exit (end of segment 
i = 8), providing the corresponding flow values are computed by counting the number of ve¬ 
hicles (irrespectively of being connected or not) that cross the virtual detector locations during the 
corresponding period. 

The estimation scheme is configured in order to estimate the mainstream densities and the total on- 
ramp flow. In the first scenario, we suppose that all vehicles are connected, thus the segment speed 
is computed using information from all vehicles present; whereas in the second scenario, variable 
percentages of connected vehicles are assumed to be present in the highway, thus the estimation 
algorithm uses a less accurate measurement of the segment speeds. In the latter case, vehicles 
entering the network are randomly marked as connected, according to the assumed penetration 
rate and based on a uniform distribution. 

The ground truth, used for evaluating the estimation results, is represented by the total density in 
each segment and the on-ramp flow. The segment densities pi{k) are computed by counting the 
number of vehicles that are present within segment i at time instant kT divided by the segment 
length (0.05 km); whereas, the on-ramp flow r^ik), is computed similarly to the computation of 
flow by a virtual detector, i.e. by counting the number of vehicles leaving the on-ramp (Lane 7 in 
Figure 3) and entering the mainstream during time interval (k, k-\-\\. 


Performance Evaluation Under 100% Penetration Rate 

In the first scenario, it is assumed that all vehicles are connected, thus capable of communicating 
their speed and position to the central authority; consequently, the filter is fed with accurate mea¬ 
surements of the average segment speeds. The parameters used for the Kalman filter are the fol¬ 
lowing: Qi^i = 1 , 1 < z < A^; Qi^i = 0.01, N <i<N%lr + ls-, Qij = 0,y i ^ j; R= 10 ^ ; 

/i = (40, ...,40);H = /(^+;^+;^)x(zv+/,-rzd- 

It should be noted that under 100% penetration rate, the filter would not be really needed, be¬ 
cause the total densities and flows could be directly calculated based on vehicle position report. 
Nevertheless, we apply the estimation scheme also for this case to enable comparison with lower 
penetration-rate results. 

The following performance index, formulated as the Coefficient of Variation (CV) of the root mean 
square error of the estimated density p with respect to the ground truth density p, is proposed for 
numerical evaluation of the estimation methodology: 


CVp = 




( 20 ) 


The obtained results are characterised by CVp = 14.9%, which imply a very accurate estimation; 
this is also visible in Figure 4, where it is clear that the produced estimates capture with high 
accuracy the dynamic variations of traffic states. The estimation of the ramp flow suffers from 
some delay that is caused by the distance of the downstream detector from the on-ramp location; 
in fact, the dynamics of flow entering from the on-ramp can be captured only after the on-ramp 
vehicles reach the measurement location (i.e., at the end of the stretch). As a matter of fact, these 
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FIGURE 4 The trajectories of real and estimated densities for the NGSIM case study, as¬ 
suming that all vehicles are connected (perfect knowledge of segment speeds). 


discrepancies may be substantially reduced by placing an additional flow detector immediately 
downstream of the merging segment. 

Performance Evaluation Under Various Penetration Rates 

Within corresponding experiments we consider mixed traffic (of connected and ordinary vehicles) 
with different penetration rates of connected vehicles. Essentially, this corresponds to evaluating 
the robustness of the approach to potentially less accurate measurements of the segment speeds. 
In this case, because of the limited number of vehicles transmitting information, two main issues 
may arise: first, the measurement retrieved from a limited number of vehicles may not be fully 
representative for the real average speed of a segment; and, second, there may be no connected 
vehicles within some segments when the measurements are retrieved. Potential reasons for the 
former are related to microscopic and behavioural phenomena; for example, consider the case 
that there is only one connected vehicle in a segment and its speed is retrieved after the driver is 
forced to strongly decelerate because of an aggressive lane-changing manoeuvre from some other 
vehicle, or the connected vehicle is driven by a particularly fast (or slow) driver. To compensate 
for potential large errors in the reported segment speeds in such situations, a moving average of the 
reported segment speeds is employed for feeding the filter, rather than using instantaneous segment 
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Penetration rates 

FIGURE 5 Comparison of performance indexes of the relative error for the NGSIM data, 
considering different penetration rates of connected vehicles, using the current speeds and a 
moving average of the reported speeds during the last 10 s. 


speeds. In particular, the last three measurements are averaged (thus, considering a time window 
of 10 s), ignoring measurements at time steps when no speed is reported (because no connected 
vehicles are present). 

The parameters used for the Kalman filter are kept the same as in the previous case. From the 
performance index comparison displayed in Figure 5, it appears that the use of a moving average 
for the reported segment speeds is indeed advantageous over the use of instantaneous speeds for 
low penetration rates and that the performance of the estimation algorithm slowly degrades with 
respect to a decrease of penetration rates, with the performance index CVp never reaching 35% in 
all tested cases. Thus, the estimation scheme is seen to produce good results even if the penetration 
rate of connected vehicles is very low (e.g. 2%). For the sake of brevity, only the results related to 
one of these tests, namely with a 5% penetration rate of connected vehicles and when feeding the 
filter with the moving average of the reported speeds during the last 10 s, is illustrated in Figure 
6; despite the obviously less accurate estimation, compared to Figure 4, we can observe that the 
estimator is capable of following the real density pattern and the on-ramp flow dynamics (with 
some delay). 

We performed also some tests in order to evaluate the magnitude of the error related to the incorrect 
measurement of average segment speeds for some cases where low penetration rate of connected 
vehicles is considered. In particular, we run 10 simulations considering different sets of vehicles 
being connected (for a fixed penetration rate), computing then the mean covariance w of the re¬ 
lated conservation equation error due to discrepancy between the segment speed computed using 
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FIGURE 6 The trajectories of real and estimated densities for the NGSIM case study, as¬ 
suming a 5% penetration rate of connected vehicles; the filter is fed with a moving average 
of the reported speeds during the last 10 s. 
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information from connected vehicles Viik) and the real average segment speed vtik), according to 

^ = ■ ( 21 ) 

’ i=lk=l 

The resulting value for a penetration rate of connected vehicles equal to 20% is w = 8.1 veh^/km^, 
while for a penetration rate of 5% the resulting mean covariance is w = 19.2veh^/km^. These 
values suggest that better performance may be achieved by utilising higher values of Q for lower 
penetration rates. Nevertheless, experimental tests not reported here) indicate that, for low pen¬ 
etration rates, the performance of the proposed estimation scheme does not vary significantly for 
different values of Q. 

Unfortunately, NGSIM data were collected for a short highway stretch within a peak period char¬ 
acterised by congested traffic conditions only; therefore the capability of the proposed estimation 
scheme in capturing the flow breakdown occurrence cannot be tested with these data. In order to 
extend the evaluation range of the proposed estimation scheme, another case study is performed 
and it is presented in the next section. 

Case Study 2: Highway A20, the Netherlands 

Network and Data Description 

The second case study is based on real data obtained from detectors of a stretch of the highway A20 
from Rotterdam to Gouda in the Netherlands, taken from {34). The topological characteristics of 
this network, which incorporates a non-trivial combination of lane-drops, on-ramps, and off-ramps, 
the congestion pattern, and the relatively closely-spaced detectors make it a challenging test-bed 
for the estimation scheme proposed in the previous section. 

The considered stretch is about 11 km in length, includes 4 on-ramps and 3 off-ramps (however, the 
on- and off-ramp “Gas station” are ignored since they are characterised by extremely low flows). 
The stretch includes a lane-drop located at 4477 m (see Figure 7). There are 32 detectors placed 
rather homogeneously at a distance of 300 m on average, as illustrated in Figure 7. In addition, a 
flow detector is located at off-ramp “Moordrecht”. 

Data from the morning rush hour of Monday, June 8, 2009 are employed, where a strong congestion 
is created around 6:20 AM because of the increased flow entering from on-ramp “Nieuwerkerk a/d 
IJssel”; consequently, the congestion spills back and worsens at the lane-drop area, reaching up to 
the “Gas station” area; then it disappears after 7:40 AM because of the reduced demand. A 3d-plot 
illustrating the traffic conditions from 5 AM to 9 AM is shown in Figure 8. This congestion pattern 
allows to test and evaluate the proposed estimator under varying traffic conditions, which include 
the formation and dissipation of a stretch-internal congestion which is not visible at the stretch 
boundaries. 

The network is space-discretised with A = 31 segments, where each segment is delimited by a pair 
of detectors. According to this space-discretisaton, on-ramps are placed within segments 3, 17, 
and 25; whereas off-ramps are placed within segments 14 and 21 (see also Figure 7). 

The proposed estimation algorithm is fed with the following information: flow measurements 
retrieved from a limited number of detectors according to the configuration shown in Figure 7, 
that is sufficient to guarantee observability as explained in the previous section (for details, see 
Bekiaris-Liberis et al. (27)); and speed measurements retrieved from all detectors, on the basis of 
the assumption that similar information may be obtained from connected vehicles reports. Specif¬ 
ically, the detector located at the highway entrance (0 m) is used to obtain the input of the system 
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FIGURE 7 A graphical representation of the considered stretch of highway A20 from Rot¬ 
terdam to Gouda, the Netherlands. Detector positions are indicated as the distance (in m) 
from the network entrance. The detectors used by the estimator for obtaining flow measure¬ 
ments within the case study are underlined. 
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FIGURE 8 The space/time evolution of real densities for Case Study 2. The presence of a 
congestion originated at the Nieuwerkerk a/d IJssel merge area that spills-back and further 
deteriorate at the lane-drop, marked as a dotted line, is clearly visible. 
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qQ\ all ramp flows are assumed unknown, therefore four additional detectors are utilised (one be¬ 
tween every pair of unmeasured ramps), choosing the ones located at 4477 m (segment i = 13), 
5125 m (segment i = 16), 6566 m (segment i = 20), and 8162 m (segment i = 24); and, finally, one 
measurement is taken at the network exit (10828 m), employed as (see also Figure 7). 

The ground truth is represented by the densities computed using the measurements from each 
detector as p = q/v, where q and v are the measured flow and speed respectively; and by the only 
ramp flow measurement available, namely the flow at off-ramp Moordrecht. 

Using the same measurements both for feeding the estimator and for constructing ground truth 
implies that the estimator is not subject to any measurement error. However, in order to assess 
the performance of the estimation scheme in more realistic scenarios, two additional cases are 
considered. In the first, a zero-mean Gaussian white additive noise, characterised by a standard 
deviation of 300 veh/h, is added to the flow measurements that are utilised by the estimator; in the 
second, a zero-mean Gaussian white noise, with a standard deviation of 5 km/h, is added to the 
speed measurements that are employed by the filter, in addition to the flow measurement noise. 

Performance Evaluation 

The parameters used for the Kalman filter are the following: Qi i = 1, 0 < / < //; = 0.01, N < 

^ ^+ir + ls'',Qi,j = 0, V/ 7^ y; /? = 100/(;^+;^)x(z,.-f4); M = (4, • • • ,4); // = /(yv+/^+4)x(iv+z,.+4)- The 
performance index in Equation 20 is used for numerical evaluation of the estimation methodology. 
In all the tested cases the estimates are very accurate. In the first case, where noise is not added to 
the measurements utilised by the estimator, a performance index value CVp = 14.6% is obtained. 
When the flow measurements used by the filter are subject to noise, we obtain CVp = 15.0%, 
which is only a slight degradation with respect to the noise-free case, since the proposed estimation 
scheme is capable of efficiently handling the effect of additive noise. On the other hand, when noise 
is also added to the speed measurements, a more significant deterioration of the performance index 
(CVp = 16.1%) and noisier estimates are obtained; this is due to the fact that, in this case, noise in 
not purely additive, but it also appears as multiplicative within the dynamic equations of the filter. 
By visual inspection of the 3d-plot shown in Figure 9, which illustrates the estimated densities 
when noise is added to both flow and speed measurements, we can see that the estimator reliably 
reconstructs the congestion pattern within the network, both in time and space. Furthermore, 
Figure 10 depicts the measured and the estimated density trajectories around the location where 
traffic congestion starts: we can see that the estimation scheme captures the onset of congestion 
with accurate timing and at the correct location. Finally, Figure 11 illustrates the results of the 
flow estimation at off-ramp Moordrecht (the only one for which flow measurements are actually 
available), from which it is clear that the algorithm estimates the off-ramp flow with good accuracy. 

CONCLUSIONS 

The macroscopic model-based traffic estimation scheme proposed in (27) has been tested on two 
case-studies using data stemming from different real experiments. In both cases, the results demon¬ 
strated the effectiveness of the methodology, both in qualitative and quantitative terms. The first 
case, based on the NGSIM microscopic data, demonstrated the capability of properly estimating 
the traffic state using aggregated information retrieved from moving vehicles; a set of experiments, 
using different penetration rates of connected vehicles, demonstrated the robustness of the pro¬ 
posed methodology with respect to potentially inaccurate speed measurements. The second case, 
where a longer highway stretch and time horizon are employed, has permitted us to verify the ca- 
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FIGURE 9 The space/time evolution of estimated densities for Case Study 2, when additive 
noise affects both flow and speed measurements. By comparison with Figure 8 it is clear that 
the estimator is capable of reproducing the real congestion pattern (both in time and space). 
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FIGURE 10 The real and estimated density trajectories at the location where congestion 
starts, for Case Study 2, when additive noise affects both flow and speed measurements. The 
proposed scheme is capable of estimating with high accuracy the onset of the congestion. 
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Off-ramp Moordrecht 



Time (h) 

FIGURE 11 The measured and estimated flows at off-ramp Moordrecht, for Case Study 2, 
when additive noise affects both flow and speed measurements. 


pability of properly estimating the traffic state in case a congestion is created within the network, 
making this estimation scheme particularly useful for traffic control applications. 

Ongoing work involves testing more complex scenarios, using microscopic simulations, to assess 
the effectiveness of the approach in case connected vehicles are characterised by a significantly 
different behaviour; e.g., when they are automated. 
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